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ABSTRACT 


Pau rning increases base pressure. The character- 
istic and linearized equations with heat addition are discuss- 
one scOnaltrtom tor maintaining irrotationality with 
heat addition 1s derived. A transition model, which extends 
euenGroceo-Ucces base flow model, and a spin model are dis- 

Bu c paed results are compared with experiment. The 


emona location for supersonic external burning 1s given. 
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NOTATION 


English 

a speed of sound, feet/second 

Cr skin friction coefficient, dimensionless 

Co ec lC hoat atus stani pressure, 
foot-pounds/slug-degree Rankine 

F dimensionless ratio of various boundary 
layer thicknesses 

h S Il lc Ont hal, 1 Oot -pounds/s lug 

k ars layer Mixine Coetiicient, dimensionless 

E Dove bench rest 

M Mach number, dimensionless 

m mass flux, slugs/second-square foot 

n Donma Coordinate, Leet; also 
C oc (y profile integer 

p static pressure, pounds/square-foot 

Q Deu added per unıt mass, foot - pounds/slug; 
Hucexdrmenssonless heat addition rate of Eqn. (12) 

q eat addition rate, foot pounds/slug-second 

R Universal gas constant, foot-pounds/slug-degree 
Rankine 

Re Reynolds number, dimensionless 

r mgl coordinate, "eet 

S unc ycwentropyj9toot-pounds/slug-degree Rankine 

S streamline coordinate, feet 

T S Ií temperature, degrees Rankine 

U velocity vector, feet/second 

U,u speed, feet/second 








w/r 


z/r 


Greek 


Q 


ES 


AS 


m 


=| @ D 


(Y 


speed along streamline, feet/second; also 
nondimensionalized by stagnation speed 
of sound, a. 


dimensionless heat zone length 


amens ironless position on body or shear layer, 
measured from base 


flow field type, dimensionless 

(M2 - 1)? 

ma io of specific heats, dimensionless 
boundary layer thickness, feet 

bpousmdgary layer displacement thickness, feet 
paumdary layer momentum thickness, feet 
dimensionless small quantity, 0 <e << 1 


Ms ri iro characteristic coordinate, feet; 
also dimensionless heating parameter of Ref. 7 


flow angle, radians 


dimensionless ratio of average internal to 
ese nal speed 


Mach angle, radians 

r ili] Meyer angle, radians 

Munning characteristic coordinate, feet 
(5 1 sluss/cubic foot 

perturbation velocity potential, square feet/second 
pu MB radıans/second 


angular velocity, radians/second 


Superscripts 


( )' dimensionless perturbation quantity 
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subscripts 


1 


er 


Ert 


average condition 

ficestream condition 

Dase condition 

condition at edge of shear layer 
laminar 

olal or stagnation condition 
turbulent 

laminar transition point 


pupbusent transition point 
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I. DNRRODUCT ION 


Projectiles in supersonic flight experience drag. One 
Eource of drag is due to the pressure acting on the base 
| (rear a e oL the po) Lile which, under normal cir- 
cumstances, is below ambient pressure. The base drag force 
CR oual to the base surface area normal to the direction 
s ene multiplied by the mean difference between the 
doent and base pressures. 

External burning has been studied by many authors [Ref. 
MIUEFoPdetermine its ability to reduce the drag of a body or 
mo produce thrust. “Ihe basic objective of external 
Durning in A Sie SOC Llon mre, TS toi create compres — 
sion waves which impinge on the base flow region and increase 
we aC pressure within the recirculation zone, thereby 
increasing the pressure acting on the base of the body to 
meduce base drag or produce thrust. 

Other methods of reducing base drag or producing thrust 
Wmemavallable. The rocket assisted projectile provides thrust 
from combustion of self-contained fuel and oxidizer. With 
Perera burning, additional fuel Can replace some of this 
Borken Since oxygen is used from the atmosphere {for burnèng 
Ehesexcess fuel. The ram jet assisted projectile delivers 
excellent performance but does require air inlets and has 
additional hardware costs. The method of base burning offers 
advantages similar to external burning, but has not proven 


mecmabrility to create thrust. 
12 
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The external burning assisted projectile provides a 50% 
inerease in range when fuel rich propellant is substituted 
for one third of the conventional projectile payload. In 
addition to an increase in range, external burning can be 
Heed to decrease the projectile flight time to the target, 
miebease 11S kinetic en@rgey at the target and provide a 
flatten trajectory. 

The base pressure ratio, Py/Po > is a measure of the ef- 
memmvencss of the extemal burning process. As shown in 
Pig, 2, the behavior and effect of this ratio can be related 
Bostene low angle, 0, of the boundary layer streamline at 
De basc. 

purposes of analysis, the exteinil burning process 
may be described thermodynamically as a heat addition term 
Me energy equation, bearing in mind that numerous other 
factors must be considered with actual combustion, such as 
@iemetteect or pressure on combustion rate. The linearized 
problem involving two dimensional planar uniform flow with 
mue cStos3ddytion has been solved by Tsienand Beilock in 
mE U Or Ref. 2,3]. In supersonic flow it was shown 
MINE Eoessure perturbatjon, due to the heat souree, is 
mud miNDed along a characteristic. Further study of the 
linearized flow equations by Rues has shown that there 1s an 
equivalence between heat, force and mass sources [Ref. 4]. 
These unit sources were given an integral formulation by 
Fuhs to determine the flow field perturbation resulting from 


mESSuroe zone [Ref. 5]. It is tlh study of these linearited 


14 
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4 


results which suggested that for experimental investigation 
the external burning process may be simulated by supersonic 
uude Wall contours which gencratc identical charactéristics. 

Heat addition in a two dimensional] planar supersonic flow 
Was investigated by Strahle, using one dimensional and lin- 
earızed heat addition theory, and shown to reduce base drag 
Piepaeoduce.thrust (Ref. 6,7]. In this study, the base flow 
va described by the Crocco-Leés |Ref. 8] model, and the 
urs base pressure ratio was determined for heat addi- 
peones of different intensity, location and extent.  Ex- 
perrmental investigation [Ref. 9-11] has demonstrated the 
feaotorlity of external burning and shown that it can in- 
Mas base pressure. Heat addition in an axisymmetric 
Supersonic flow has been studied by Hosack and O'Leary using 
n c (hod of Characteristics and thé Korst base flow model 
Mew. 11]. 

M object of thé present investigation was to study 
base pressure in two dimensional planar and axisymmetric 
supersonic flow with external heat addition using the method 
Dna xgctorrstics and the Crocco-Lees base flow model. The 
cffects of varying Reynolds number, including boundary layer 
mansion. and projectile spin havc been included. A con- 
current experimental investigation studied the effect on 
base pressure of various axisymmetric nozzle contours which 
Semllated he external burning process. 

It should be mentioned that several combinations of ex- 


ternal heat addition analysis and base flow models have not 
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been investigated theoretically. These include two dimen- 
sional planar method of characteristics and linearized theory 
mech cine Chapman-Korst base flow model, and linearized axi- 
symmetric flow with both the Chapman-Korst and Crocco-Lees 
base flow models. The linearized axisymmetric flow problem 


Mimi be discussed later. 
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EE 


The theory employed for this investigation must describe 
NI USseady flow in two entirely different flow regimes, 
Emi —hecexternal flow is supersonic and can be accurately 
ucc DY the method of characteristics. The internal vis- 
EN NMIoc OF recirculation zone, is dissipative and diffi- 
muB to describe. As in the two dimensional case, it was 
ee amed that the Crocco-Lees model [Ref. 1,12-15] adequately 
meeeribes the axisymmetric internal flow. Boundary layer 
mapes (laminar or turbulent), thickness and transition were 
Exeuded as a function of Reynolds numer since these effects 
Mmeoroundly influence the base flow behavior. 

Meme maracteristic equations have been derived for flow 
wheat addition (Ref. 16,17]. In vector notation the 
wady equations for flow without viscosity, heat conduction 


Enebody forces are 


UU 0 (1) 
pU-VU + Vp = 0 (2) 
WEE dU) sq (3) 
De zo! (4) 


Erere p 1s density, U is velocity, p is static pressure, h is 
eut@malpy, q is the time rate of heat addition per unit mass, 
e ij universal gas constant and T is static temperature. 
Ma natural coordinate system [Ref. 18], Fig. 3, the 


> 


characteristic equations were derived from the above equations 
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in the following convenicnt form (Appendix A) 








ec 9h 
o(v-8 0 b : C 05 0 
)(v-0) 5 ) = ( "Hm. - Ver ) s Till ss Sar (TS 5 urn) (5) 
p 
= nu dh 
nero) 2 asind _ b š cosy as | 0 
dE SCH r/r, Sn ) sinu + (152 5n) (6) 


were 0, VY and y are respectively the flow, Prandtl-Meyer 
w ach angeles, U is speed, S is specific entropy and ho 


is total enthalpy, hy = h + u she left and right running 
eharacteristics are n and £ respectively. The natural co- 
Ordinates, s and n, are along and normal to the streamline. 
as two dimensional planar flow, a is 0, and a is 1 for 
ERusumetrijc flow. The characteristic equations, (5) and 
K have becn multiplied by the characteristic length of 
base radius, rj, to obtain a dimensionless form. Therefore, 
in contrast to the dimensional coordinates n, & and n used 
EN D 1 c:neseoordinates n, E and n in Eqns. (5) and (6) 
Me understood to be dimensionless. Also, the specific heat 
s Constant pressure, EE was assumed constant in the enthalpy 
relation h - Ge 

IN eaten is tie computetions w1th Eans. (5) and (6) 
were simplified by neglecting the term pa : un )/T, , 
which is the vorticity. One may assume the external flow to 
be irrotational before heat addition and then determine the 
RER on heat addition which preserve irrotationality. 


Meco sSeuncorem for steady flow may be written (Appendix B) 


TVS + UXQ = Vh, (7) 


Mere 2 = VXU. Taking the curl of Ran. (7) 
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— 


Uevn = 2-VU - QV-U * VTXVS. (8) 


Assume 2 = 0 before and during heat addition. Then, from 


Eu 08). the condition which preserves irrotationality is 


VTXVS = 0. (9) 
fome redefines the heat addition and has uniform upstream 
Sonditions at infinity [Ref. 19], then 


V(hy - Q) = 0 (10) 


where Q is the heat per unit mass which has been supplied to 
BH fluid element. Using Eqn. (7), with 2 = 0, Eqns. (9) 


ama) (10) yield the following result 


VOXV(4U*) = 0 (Ti) 
E ors the condition preserving irrotationality; that is, 
the gradients of kinetic energy and heat TIE must be par- 
EN t O-ys-a small perturbation, (Q/h,) = €, then Eqn. 
Meets of order e and the irrotationality condition is 
EE EE Satistied to order e. To neglect the vorticity 
mom in Eqns. (5) and (6), one must insure 

NE Uv 


T i "bUc T poc SMS (12) 


XO 





Tmo weommumattols discussed in Section V, the quantity 
i bon. (12), oe ee was assigned a constant value 
independent of local pressure, within the heat zone and was 
Men as zero outside the heat zone. Thus, within the heat 
Es /0as- approximately parallel to the gradient of kinetic 


eieregyeand the irrotationality condition of Eqn. (11) was 


Zi 








Euustied. This condition was not satisfied along the inner 
Rncuter streamlines of the heat zone, since a jump in Q 
weed normal to the streamlines. However, by assigning a 
Matue for Eqn. (12), (ry) / (Uc. T) ge l] the irrota- 
meondlity condition was satisfied to order epsilon along the 
a C imline boundaries of the heat zone. 


EquUataons (1) through (4) may be combined to give 
ues — 0 O) = (13) 


Lucre a is the speed of sound. Assuming small perturbation, 


mon. (13) may be linearized (Appendix G) 





29 - m 28 -g 14 
$ = ox? - q ( ) 
meee rm 1S the Mech number at infinity, q' is the heat ad- 


meno, q — h q', $ 1s the perturbation velocity potential, 
NEU” + 0", UY = V$ and T » T. (1*T'). In Cartesian co- 


CO 


Exwnates, Eqn. (14) may be written 


ER wu a (15) 


where 8* = MZ - 1. For supersonic flow, the solution of 
Eqn. (15) is found by integrating the sources in the fore- 


cone volume, Ve Mer. 20], 
ue f Sc ss 16 
a Y ERAN ON RN P 
f 


opa stne linearized theory was not pursued further, it 
piamua be noted that this theory could be used instead of 


a cethod of characteristics to establish the effect of 
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Eunal burning at the base flow shear layer. Studies of 
Finearized supersonic flow through ducts have provided solu- 
Mons for an annular source distribution [Ref. 21-23], which 
mud be used in conjunction with centerline sources and 
Ems [Ref. 18] to describe the axisymmetric external flow 
with burning. 

lbs tuo djimenssonal planar Crocco-Lees theory [Ref. 8] 
considers the base flow as a separated boundary layer merg- 
Ea uto a free jet profrle and derives integral equations 
memcescribe the mean values of the flow quantities (Appendix 
Brom the integrals of the conservation equations, a 
le ordinary differential equation is derived. This or- 
Misery differential equation has a saddle point singularity; 
ss ce of “the singularity is a key feature of the Crocco- 
Lees model. The behavior of this equation is partly deter- 
Ned by an empirical F - K relation and by the nature of 
Menear layer mixing. k is the ratio of average internal 
EL ow Velocity to the external velocity, u,/u,. Fis a 
metio formed by using the various boundary layer thicknesses. 


Mm turbulent shear layer mixing coefficient, k was as- 


de 
wehe Constant value, 0.03, by Crocco-Lees. The laminar 
Male used for this investigation was 0.0033, about 1/9 the 
merbulent value. 

conceda ta profile form fae@tor; when the*veloc- 
profile changes shape, the value of F changes. As dis- 


u ein reference 8, the PF - k functional relationship 1s 


based on many known viscous flows. For example the Blasius 
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Bon 1S represented by a point in the E k plane. 
meee r-oKkKan boundary layer solutions, which involve either 
emor able or adverse pressure gradients, become a curve in 
EN D Ke plane, The velocity profiles of jets under cer- 
conditions are described by the error function; these 
MES become a point in the F - k plane. Likewise the Howarth 
su Quy layer solutions for flows in pressure gradients are 
Per ve in the F - k plane. 

he” PF - K curve corresponds to viscous flow near walls, 
E wboeundasy layers, as well as free jets. Furthermore, 
memerse, zero or favorable pressure gradients are represented 
BEN s:Uutrons on the F - k curve. Heat addition in a super- 
sonic flow generates compression waves. When z/r < 0, the 
compression Tv eCSmumpance waincad sot the base; the boundary 
meer om the body and the shear layer experience an adverse 
Pressure gradient. When z/r > 0, the compression waves im- 
pinge downstream of the base and only the shear layer flows 
EUN be adverse pressure gradient established by the exter- 
EN "uu. 51znce the P - k curve incorporates pressure 
E ents; the presence of an adverse pressure gradient due 
Eueatcadditron should not alter its validity. As stated 
ENrclcrence 8, the F - s Meus additional experimental 
Location; this applies to the case with external burning. 
er discussion of the F - k curve may be found in refer- 
wes 6, I2, 13, 14 and 15. 

Mea econ mulation of the Grocco=Lees theory for axisym- 


eow [Ref. 1,12-14] resulted an the same single ordinary 
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differential equation which had been derived for planar flow. 
In two dimensional planar flow this equation is numerically 
[miteeerated by starting at the singularity, or so called 

u cal point, and working toward the base, using the 

E mndtl-Mceyer relation to obtain flow angle for various Mach 
Ambers. With external heat addition, one may no longer use 
EL xandtl-Meyer melation for flow angle at the shear layer, 
mimiec the heat addition, described by either linearized or 
ExrEduimensional] flow, introduces a correction term [Ref. 6,7]. 

NX uxssedusescannot be used for the axisymmetric prob- 
eit hout heat addition due to the 1/r term arising in the 
NUnüusty equation; thus, one must start at the base and 
mimeeerate toward the critical point. As was anticipated, 
meeempts Co numerically integrate through this saddle point 
memeularity were unsuccessful; however, by repeated integra- 
tion one may determine two Mach numbers on the downstream 
Side of the expansion fam which straddle the saddle point 
EuxobLarity, thus establishing the base pressure within the 
merits given by the two Mach numbers, Fig. 4. This technique 
D cused with Egns. (5) and (6) for both axisymmetric and 
two dimensional planar flow. 

Pominecorpemavwe tlie effects of various boundary layers on 
miemoody (laminar, transitional, turbulent) and shear layer 
EUM tIono rt was necessary to formulate a transition model, 
n >. Whiihe this model treats the boundary layer mass flux 
Sete rentiy, it is not a departure from the existing Crocco- 
E heory. In fact, this transition model extends the 


moco Lees base flow model. | 
Du 
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“mce the transition from laminar to turbulent flow is 
more or less probabilistic in nature and, depending on the 
Reynolds number, may occur on the body or along the shear 
Mier, it wis assumed that the transition process could be 
@eceribed by an integrated normal distribution and two 
transition Reynolds numbers. The transition Reynolds num- 
Ee Eat ke dt the start.and.end.of the transition process, 
mathe integrated normal distribution ranged from 0.02 to 
uon the span between the two specified transition 
Reynolds numbers. This assumption replaces the previously 
assumed transition curves [Ref. 8]. 

li estransrtron model must provide several parameters 
Nace required in the Crocco-Lees base flow formulation. 
am; sl establish the values for the mixing coefficient, k, 
along the shear layer and the boundary layer mass flux at 
ma base to permit determination of F. 

The boundary layer was assumed to be at constant static 
E ssure and total enthalpy, and the velocity profile was 
assumed to have a 1/n dependence, where n is an integer, 

s - 7. This lormulation may be integrated (Appendix D) 
fOrspoth planar and axisymmetric flows to give the boundary 
E sssoflux at the base as an explicit function of veloc- 
mepeprorile, i/n, base radius, Yu» and boundary layer thick- 
D if the flow over the body is completely laminar, 
spa ven by the compressible Blasius solution [Ref. 8], 

Sy = (5.041. 2(y-1)M?)L/(Re,) *, ansehe value of n = 1 is used 


ue velocity profibe [Ref. 24]. If the boundary layer is 
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eommbetely turbulent, the n is taken as 7 and 6 is computed 
Prone the compressible one-seventh power solution [Ref 8], 
E - (0.037,/8ey) ^!) (6g /Cg ) 6/699). The term (Cgy/C¿.) 
ena ratio of skin friction E icuts LOr compressible 
wM ompressible turbulent flow [Ref. 24]. It was recog- 
Meana the boundary layer thickness computed for planar 
flow will be different from the axisymmetric thickness. 
This difference will be small if 6/7) and radial velocity, 
wn are of order e. The dimensional term pv/r from the 
axisymmetric continuity equation will be of order e* upon 
nondimensionalizing with ó and U. Thus, the governing equa- 
Bons anda their solutions are identical to order € in both 
eases. 

The first transition Reynolds number, Re,» is where the 


Eu :7tlowmstarts transition to turbulent. The second 


transition Reynolds number, Re, po is where the flow must be 
eompletely turbulent (Reir : 4Re, ,.)- AS one changes values 
of base Reynolds number from Re, to Re, a> the values for n 


End 8, which determine the boundary layer mass flux at the 
meee, Must vary from 1 to 7 and 91 to T respectively. This 
W a on was computed by the error function, erf, since it 
is readily available in FORTRAN and can be directly related 


mee thewintegral of the normal distribution. A numerical 


Enos cor 4s which takes the value of erf from -0.96 to 


IM O qS used to normalize the difference Cer - Re, ..). 
E SR = D 
Det A 4 (Re... + oe Then 
N 5 
` 4 = 1. oil m D = = d 
ERE, pré. Mc ROSA a pen (17) 
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So. when Re, > Re, y 


S ASE - 9, C: t ERF,) + Sp Cs - ERF y) (18) 


and 


e lee ER.) € 75 - ERF,) (19) 


were the values used to compute the boundary layer mass flux 
meee base (Appendix F). 

When transition from laminar to turbulent flow occurred 
Mene ce i layer, the value for the mixing coefficient, 
atra particular point was computed in a similar fashion 
Me the error function to transition from laminar mixing, 
k 


Imanu rbiulentrinixinge, k linus. 


I^ T: 


M eor e hear layer mixing om base pressure is 
meun schematically in the momentum control volume of Fig. 7. 
iiewaxial pressure on surfaces cd and de, and momentum flux 
lane the volume, de, must be balanced by the pressure on 
surfaces ab and bc and momentum flux entering the volune, 
K thus, £0 Maintain this balance, an increase of momen- 
ALO: along cd due to mixing requires a corresponding de- 
c se or base pressure, Pp > along ab. 

One additional assumption was necessary to obtain good 
Ewcementowrth experimental data. At low Reynolds numbers 
and completely laminar flow, it was necessary to increase 
micevalue of the laminar mixing coefficient, ki, in direct 
u I I On O the :rncrease ot the thickness of the boundary 


Ext oche base. This assumption may be expressed as 


E 
k, 7 k(Re,, /Re, ) 7 (20) 
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where Re, is the transition Reynolds number, Re, is the 
eos number of the body, based on length or chord, L, 
and k is the laminar mixing coefficient when Re,r = Re, 
s justification may be provided for Eqn. (20) by consider- 
mmeeche Blasius boundary layer solution where the mass flux, 
m, is given by 3.28pux/ (Re) and the Crocco-Lees theory 
where kpu is taken equal to dm/dx. If one forms the deri- 
vative of 3.28pux/ (Re) 5, dm/dx is found to vary as 1/ (Re), 
BEcby suggesting the form assumed in Eqn.: (20). Although 
k... Berg come +ant at 0.03 for all calculations, a similar 
Ames is using the incompressible turbulent solution suggests 
that Kr might vary as 1/ (Re) V5. 

Eu 1 -osuspespedP taut tbe effect ot projectile spin 
REuPO be included 3n the mixing coefficient [Ref. 15], k, 


By using k = k(1+rw/U), where ru/U is the spin number. 


spin 


However, since k always increased with an increase in 


spin 
Spin number, and, since it has been determined through numer- 
Ee experiments With the Crocco-Lees theory, Fig. 8, that 
increasing k decreases the base pressure, this simple hypoth- 
ees) CONtradicts experiment [Ref. 25,26] where, in some 
cases, increasing the spin number increased the base pressure. 
ie. 2 tos clear that at certain values of Reynolds num- 
Der, Rey, an increase of Reynolds number resulted in an in- 
Base in base pressure. This fact suggests that the spin 
number could be used to shift the pp versus Re, curvo, 
Emu) oaxnstead of increasing the mixing coefficient. Sev- 


Ecefthods were attempted at this point, such as increasing 


Crt 
Y? 










Reference: k = 0.03 
ern? 
EU n-/ 
Mach = 2.0 
r r rn on e 
U.O = ——— 2/8 
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D. u minine coefficient in 
OO shear layer ME 
NI = Fatıo of base radius to 
boundary layer thickness 
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n = exponent for boundary layer 
velocity profile 
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Figure 8. Sensitivity of Base Pressure Ratio 
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Re, By [Ref. 27] (1*1) ) * or (4€ C3 *) , but the method 
which provided the best correlation was to simply decrease 
the transition Reynolds numbers. The spin transition 
Reynolds numbers used to compute Py / P versus Re, were form- 
ed by dividing the corresponding no spin transition Reynolds 


numbers by (Or Gol), RE of base pressure data 


£ Ca spinniner projectile 1s .given in Section V. 





Pele te, PeRIMENTAL INVESTIGATION 


momearry cut experiments related to the five inch pro- 
Tooele, the mid point of a maximum range trajectory was 
s (cq. This flight condition established an altitude 
DES 000 feet and a Mach number of 2.0 for determining 
Mita ty parameters, For the five inch projectile, this 


6 


gave a body Reynolds number of 15.7x10° , a radius Reynolds 


number of 5 


E bulent boundary layer thick- 
K o radius, ó/r, of 0.1513. The spin at 23,000 feet 
Ln csssumed fo be the same as the muzzle spin, 250 revolu- 
NOUS per second, which gave a spin number, rw/U, of 0.1595. 
Using these nomial values, it was found that all similarity 
Mame ters could be matched with a free jet Mach 2.0 coaxial 
Mate exhausting at room pressure if the center cylinder 
NEE 0.641 inches, ô of 0.084 inches, and a spin 
w revolutions per second. As stated in the introduc- 
mean, the actual combustion process was simulated with 
merous nozzle contours which duplicated the compression 
waves generated by the external burning. 

ISS Illa trono process 1s depicted 1n Fig. 10 where the 
camine deflected by the presence of a heat zone can be 
gu ueed by a nozzle wall to generate an identical flow 
tates 

Hu" yet wc perimental facility, Figs. 11 and 12, was 
EOwwmtronal in most respects and has been previously re- 


Borvegd [Ref. 28,29]. The use of cast aluminum filled epoxy 
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(b) Identical Flow Field With Compression Generated by 
Nozzle Wall to Simulate Heat Zone 


Figure 10. Heat Zone Simulation With Nozzle Wall 
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Figure 11. Photograph of Free Jet Experimental Fac 
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‘Figure 12, Schematic Diagram of Fxee Jet 
Experimental Facility 
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nozzles, Fig. 13, was the only new hardware innovation. The 
technique of casting employed a wax plug contoured to the 
desired shape with an aluminum template. This plug was 
memecd Slightiveummeaersaze and then restored to sizce using a 
amn C oat of amseneeıy e body plastic After melting the wax 
Monwi Cast nozzle, the plastic peejed away from the epoxy 
to leave a smooth nozzle wall. Without this plastic, the 
heat generated by the cpoxy had softened the wax and caused 
“ome Emer nozzle wall. 

Temo cation or the zone Of combustion to be simulated 
by the nozzle was estimated using an empirical formula for 
Be tralsonsdepreuseorsar jet normal to a supersonic frec- 
Stream weer. 350]. For Mach 2.0 flow, internal pressures 
from 300 to 500 psia were computed to provide penetration 
bo 25 inches with various orifice diameters. Thus it was 
assumed that external burning at 5 inches from the projec- 
te Centerline would be practical. Assuming a fuel rich 
propellant exhaust, ARC-168, with an available heating value 
of 7200 BIU/1b and 0.6 1b/sec flow rate, it was found using 
one dimensional theory that the flow should be deflected 
aaa 45° for a distance of 15 inches, or three diameters, 
downstream. This computation used the integral methods of 
Gëfter 5 win pave a perturbation velocity normal to the 
ono TCS Seeummer sccond. The resulting flow angle 
tangent was -163/2058, which gave a computed flow angle of 


Ec ber inner streamline of9?the heat zone. 
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Insert Showing Spacers. 


Coaxial Epoxy Nozzle. 


Figure 13 
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These estimated values for the heat zone were taken as 
imas o limat and were. used to form a test matrix and for 
um -cauentenuoz2lesdes)pno Seyeral different nozzle contours 
were Computed by the method of characteristics for 4%° stream- 
lum x ofticctron at a radius equal to a cylinder diameter 
igomstnescenterline and with different assumed lengths for 
GhcmeoubUsbmMon-zone. wBiwewbestwmnatrix:included using spacers, 
DE rt."vwi1ith eseNEnvzzle so that the first full com- 
pres ves Impingedeanead Of, at, or behind the base of 
er. FErzessuressnere measured along the base and 
EDEN rtc ane cylinder, Fig. 13, using 
Ma ON OO rS Mp eS to carry out experiments with 
Ide: were Unsuecessiul due to faulty air bear- 
INNEN H. Time daid mot permit correction of 
ts AA EU 

Segen the flow field with the stationary 
(C I S s I oz2 e C nB O seen in Fig. 14. The 
mimosa pimo tic Uniform Elow nozzle in Fig. 14a used a 
stal knife edge and the expansion fan at the base ap- 
Ia ark. The base pressure ratio was 0.635. In Fig. 14b 
the vertical knife edge was positioned to make a positive 
sss ora iL from the centerline to the right appear 
nan Neon compression, Fig. 14c, had an ex- 
pansion fan at the base but achieved a higher base pressure 
Piero Ole m to compression upstreamwef the base as 
Welieaa alone the shear layer. The compression upstream of 


ER WEE the pressure right before the expansion 
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Figure l4a, 
Uniform Flow, 
Horizontal 
Knife Edge 


| Figure 14b. 

| Uniform Flow, 
Vertical 
Knife Edge 





Figure 14c, 
Compression Flow, 


8 e | N Horizontal 
Ë BS b ` Knife Edge 


Figure 14, Schlieren Photographs at Mach 2.0 Showing 
Uniform Flow and Flow With Compression 
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maneso that 1t was greater than p_, thus the flow was near 
ambient pressure after-passing through the expansion fan. 
omo test results using the stationary cylinder are given 
able [t Ecepte tor the nozzVPe with w/r = 4.34, these 
tests indicate that compression applicd upstream of the 
Ma 272025, îs less efficient ak increasing base pres- 
sure than compression applied at the base, or slightly down- 
Stream. The nozzle with w/r = 4.34 had a strong shock wave 
ae the exit, due to compression waves reflected from the 
Ander [Ret. 15], which, with z/r = -0.5, impinged on the 
K j low region aheadioi the critical point. In turn, in- 
creased pressure across the shock wave was impressed on the 
EMbsonic flow in the recirculating zone and gave a higher 
ase Pressure, The shock wave observed with this experimen- 
enoe le would not be expected with uniform external 
Burning. 

le une experimental base pressure ratios, Py /Poo > 
w enon aS the Upper number in each box for the various 
mest positions, 2/1, of each nozzle, w/r. The cylinder 
miomucmiser, thesaxtal location where the first full strength 
Eonpsesssmon wave impinges on the body or shear layer is z 
ELE X:hcsbaseandE is the Length of the heat zone. The 
eorresponding Reynolds number, Rey > is shown below each value 
of z/r. The lower number shown in each box is the computed 


specific impulse for the assumed projectile Filo t” C Opa it roli 


DE ihe unrform flow nozzle, with w/r AS 


Macias tie reference pressure ratio to compute the 
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First Full Strength 


Compression 
hz -> 


Ps ee S -— — -— - — - => 
0 +Z 


Sketch to Table I 


Er seme | ws | wo | 05 | mo 


Re, X 10 — 95) 780 2295 SO 


Sege 1.078 LAS, Dessen 1.129 
137 169 167 158 
4. 1.20 1,14 1,085 
156 143 131 


Table I. Base Pressure Ratio and Specific Impulse 
With Simulated Heat. Addition 


lio 0.96 12055 1.074 
216 283 302 





base me suresrorcescreategd by the other nozzles. The fuel 
flow rate was taken proportional to (w/r)/(6.0), since 6.0 
was the combustion zone length assumed for the fuel flow 

99 0.6 lb/sec. 

Once Mach number and operating pressure and temperature 
wee elected, the remaining variable controlling Reynolds 
number .issthe length or scale of the apparatus. The free 
m e ey o Figure tinas sufficient mass flow to per- 
mec duplication of Reynolds number at 23,000 feet altitude 
pucvudedothecevlirnder has sufficient length. Consequently 


6 For 


the nominal Reynolds number was chosen as 15.7X10 
IO rece match Capability of the facility. 

When the design of the spinning cylinder was initiated 
membecamce apparent that a long cylinder had problems in 
care to shart whip or the first vibrational mode in 
bending. A shorter cylinder reduced the Reynolds numbers 
momeme Values shown in Table I. Conditions at a higher 
ale tude were Simulated, or else conditions for a short 
Peaeceme at. 25,000 teet were duplicated. 

ls mcsumdernseceeuldebave been) designed. Each cylinder 
Ld nsEesqguired a separate family of nozzles. A choice 
Dmsdestosuseva short cylinder for both stationary and 
SPinning tests. 

Nation nar as. shown in the sketch with Table I 
s co O olds umber slightly. The variation in Re, in 


lu s nT 1S due to this fact. 
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IV. NUMERICAL INVESTIGATION 


her computer program for this numerical investigation was 
LEustcen to calculate eather two-dimensional planar or axisym- 
Meme problems. the program used the method of characteris- 
Miles where the axisymmetric term in Eqns. (5) and (6), 
ry (sinusin90)/r, was deleted for two-dimensional planar flow. 
a IDcetrecroor radrus In the mass flux calculation, 
at the base and along the shear layer, was deleted or re- 
RENE “ine characteristic mesh, Fig. 15, was 
Rer trom the body or shear layer upstream through the 
meat zone, 3f present, to a reference freestream Mach line 
BD Xb3wwards alone a right running £ characteristic). 
Mecerror from the reference condition was used to correct vV 
imine body, since 8 was given, or 8 along the shear layer 
pence Y was computed from the Crocco-Lees theory. Each body 
em shear layer point was iterated until the error at the up- 
Stream boundary condition was less than A radians. VAt 
een port im the characteristic mesh, the value of u was 
su ermmed- tor the resulting value of v by the method given 
MA dio a lie retrospect, this method of characteristics 
conp ota tion could be made more efficient if one could devise 
a scheme which works with arbitrary points along the body or 
shear layer and increments the characteristic mesh along the 
K b O treestream Mach line. Then only one calculation 


Molcepomicdiiired along Gach characteristic, rather that the 
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iterative scheme used in the present method. This approach 
was considered, but it was not clear how one would formulate 
mre calculation for the corner or the shear layer mass in- 
fix {Appendix E). 

The flow field was calculated along the body to the base 
Emwhrchopomrnt all computed values along the last character- 
Meme were stored. These values were used to reset the flow 
or reach iteration of the Straddie method. 

Ine Crocco-Lees theory was formulated in terms of Dis 
instead of Mach number, where Be is the velocity at the edge 
of the shear layer divided by the freestream stagnation speed 
pr sound, un U/a,. To start the straddle method a value 
of w. was SIS cone downstream Side of the expansion 
man which was too large. The expansion fan was divided into 
RE EE downstream of the fan the computation proceeded 
Boward the Critical point. Each new point along the shear 
ae rewecomputed by the Crocco Lees theory (Appendixes C, E 
Endo) required several iterations along the characteristic 
EN u"cxpstream Mach Tine boundary condition, where the flow 
angle error after each iteration was used to correct the 
shear layer flow angle. 

In C en Dr tme = Wa Diane Omloved the first path 
own am Fag, 16, It can be seen that this path falis to 
ere Tight of the path passing through the Saddle point 
Seigeularity, X, and, as the value of F was decreased, the 
stope, dF/dw,, became negative at 1'. This occurred when 
the denominator used to form dr/dw, became positive [Ref. 8], 
Ic cbcoHseMe numerator and denominator of dF/dw, were 


50 








aus t) 
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Figure 16. Integration in the F - Wo 
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negative OMS transition from positive to 
negative slope with positive denominator occurred when the 
Eutial guess for w. was too Farge. By systematically de- 
@eeesing thé initial guess for Wo» Say by tenths, one 
memcually has the integration path from 2 to 2', which is 
Fani lar to the previous case, except that the negative 
elope at Z‘ 1S associated with the numcrator, used to form 
dF/dw, , becoming positive. Thus, the last two initial 
guesses for We have straddled the saddle point singularity. 
l""terval between these two values may be divided, say 
ENEUundredths. and the process repeatcd to find two values 
on Gs E EE straddle the saddle point Singularity. 
mec onc was specified to start and stop at certain 
meres of axial Position, z, and to lie betwcen an inner and 
Encre streamline. The quantity fpa/ Uc, T Was assigned a con- 
stant numerical value, between 0.15 and 0.175, within the 


o 


ea nnch gave the desired -434; deflection to the inner 
Streamline. External to the hcat zone rb q/Uc T Waste ro. 
erore deflection, the inner and outer streamlines have 
values of approximately 2.0 and 3.0, maximum body radius, 
mespeetively, <Ihesc valucs were chosen to match the rate of 
SeNpressi Ome increase used in the nozzle design. As the com- 
putatitonm proceeded downstream, the valucs of radius and flow 
pele for the inner and outer streamlines were computed to 
n cats tie boundaries for the heat zone. 


When a sufficient amount of properly located heat was 


Ep csuoCecternal flow, the basc pressure ratio was 


SZ 





computed to be greater than 1. It is interesting to note 
Mat or this condition the method of characteristics rou- 
Bine generated a compression fan at the base, rather than 
ee usual expansion fan. The compression fan calculation 
Mathematically folds the flow field on itself [Ref. 31]. 
Mats Of Course 1s physically impossible and is replaced in 
meal flow by an oblique shock. For values of lee = 1.5 
Nach 2.0, the Noss of total pressure due to the oblique 
EN cius less than 1$, whach corresponds to an entropy in- 
w ase, AS/R, of 0.01. This resulted in an error of less 
Ban 0.18" in the turning au mcs uMSUECuascneelected. There 
is also the possibility of boundary layer separation up- 
EB cam Of the base for values of pi / P greater than one 
PRet. 7]. This effect was not considered in the calcula- 


mons. The computer program is discussed in Appendix J. 





C CUS TONT TOF RESULTS 


The computed values of Py / Poo versus Re,, shown in Fig. 
9, were based on the apparent transition Reynolds numbers 


lucu the Unitorm flow nozzle experiments. These were 


: and Re, = 7.8X10°. E O 


Re & 32295210 t 


Er 
the same general features and approximate values found 
Bene periment [Rer. 32-35). It should be noted in Fig. 9 
suat for both the two dimensional planar and axisymmetric 
mow the maximum value of Di / P occurred below D and 


the minimum value of DAD occurred below Re, 


E^ 
Du dnce difterent choice of transition Reynolds 
6 6 ! 
1 re = CY = D acl e 
numbers, RO 4 y EIo ana Re, Lt pe WEM resulted in 


mnie curve of D Ya US Reynolds number shown in Fig. 17. 

One may compare the curves of Py, Poo re. Sand Pres. 17 

mele time experimental curves shown in Fig. 18. While the 
Marco means identical, they do agree in having 

the same general shape. Figure 19 shows additional experi- 
a OCO pare With Fig. 9. The locations of maximum and 
minimum computed pi, / Pa muomeonuroMbedg by the choice of values 


for Re... and Re uno itio of- for Re, ,./Re used in 


up" ETE 


these computations was based on information in reference 8. 
This ratio is an approximation which varies considerably for 
iertcrent bodies, depending on their shape and surface con- 
dition. The maximum value of Di, / Pa, was controlled by the 


Choice of ku, white the choice of Ko controlled the sécond 
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maximum value of ny De to the right of Re The slight de- 


Ert 
crease of DD to the right of this second maximum is in 
Spncementowdthethewr/ó sensitivity curve of Fig. 8, since 
imeneasing Re decreases 6 and thercfore increases r/6. 

Meee onal Curves are Shown in Fig.%@ 7; thesë are 
wansi tron and Critical point locations. The laminar and 
tumbulent tmansition points are found simply by dividing 
base Reynolds number into Re... and RO it: Tier trcal 
point location was taken from calculation with the Crocco- 
Fecomcaquaunons along the shear layer where dF /dw changed 
Menes it is the singularity 1n Fig. 16 which was discussed 
I C IOo Tn going from small to large values of 
Reynolds number, the shear layer between the critical point 
and the base gradually changes from completely laminar to 
completely turbulent as the laminar and turbulent transition 
pomme e Upstream to the base and then onto the body. 
Mere tnistoccurse. thesshear Layer maxing coefficient, k, 
gradually changes from k; ËO k... Increasing from some small 
eno Re molas number, the thickness of the laminar 
way Mayer at the base decreases until the laminar 
RAS ran pOint teaches the base and moves onto the body. 
As the boundary layer at the base transitions to turbulent 
ness increases ho thickness continues to increase 
deu bulenmt transition point moves forward along the 
Doc aa ximumiwhen the transition point is near 
Mies cane body,  Beyord@ thus point, the thickness 


decreases. 
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Inre. the value of Py / Poo starts to decrease slight- 
ly before reaching Re... amic laminar transition point 
EE on he Shear layer and the basc boundary layer is 
Emar aN hou hehe boundary layer thickness is still 
decreasing at this point and in turn should provide an in- 
Ereasc in EE the value of k along part of the shear 
Payer mi Dekween the taminar transition point and the critical 


Ponne tarts to increase from k, to Kp and consequently 


L 
dominates the overall change in pM with Reynolds number. 
The mechanism responsible for the abrupt minimum value 
of PU / P. Jiro wasemec Clear, However, this minimum 
Pocweaqleulated and 1S shown in Fig. 17 as a dashed line. 
ocio sensitivity curve in Fig. 8, 1t appeared that k 
uc becmostwiekely9"quantrty responszbTe. In numerical 
experiments this abrupt minimum was calculated by assuming 
additional mixing. The additional mixing was assumed to 
Wea normal distribution about Re 


60/2 


li EE to about ten percent of its maximum 


pes k = kp (1+ 


1. 
A EE x, was adjusted so that the 


Va luc ın 2x109 Enciso oro This additional mixing 


(rí 
was not used in any other calculations since it was not 
understood and was not clearly observed in the expcrimental 
dee uremertedrlateyrssin this section. 

Calculations were made for external heat addition in 
INoudinienctenaleplatar flow to compare with Strahle's 
pies jou Door cor cC Calculations [Ref. 7]. The values of 
Re... = 1.95x10° and Rea nt = 7.8x10° were assumed and a 
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Pave tne same value of 


Reynolds number of Re, = 20-0X10 
6/r used by Strahle. The Mach number was M = 2.3. The heat 
Ada tion zome was started at the same point, z/r = 0.4, and 
mescaened the same length, w/r = 3.0, but to avoid subsonic 
meow witch the method of characteristics it was necessary to 
increase the height of the heat zone from the value of 0.2r, 
uc x 5trahle to a value of one base radius. The values 
of Q = pp exe Lorca e accordingly to provide the 

Paes toudl Meat addition based on the rectangular heat zone 
wcd by Strahle, i.e. Q = 0.2n/tany. Ino results in lable 
Meeshiow good agreement at the low intensity of heat addi- 
Bn uuu carr agreement at the higher intensities with 
Biemietnod of characteristics program showing higher values 
one Py / Po Although the non linearity included in the method 
characteristics calculation may be responsible for part 
Pete wiiener base pressuresratio, the main cause was at- 

m ui to expansion of the heat zone. This effect was 

c lie In the bounding streamline calculation for the heat 
Zones ine tine method of characteristics program and resulted 
EuCascreater amount of total heat added than in the linear- 


zed theory calculations. 





Pineat ized Method of 
Theory l’! Na teristics 
n Q EE 
0.0 0.0 0.410 
0.4 0.167 0.669 
0.8 0.334 AS 
m 0 5 ee OZ 
Iu ih Comparison of Calculated Base Pressure Ratio with 


Eno mule dartion, Iwo Dimensional Planar Flow and Mach 
Number of 2.3. 
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Both experimental and calculated results are shown in 
Eu 04for a Mach number of 2.0. The value of P, / P - 
026314 at Re, = 10.4 was the average of nine runs with the 
Unstorm flow nozzle. The probable error about this average 
Nas t 100.004, with a maximum observed deviation of - 0.0104. 
Ncomparing the experimental values of Fig. 20 and Table I 
Annea y test results using the free jet, variations as 
NL as 7% in Pi / D Were found tor similar runs. These 
Variations were attributed to problems, such as mechanical] 
ENUumentegdurrng nozzle installation or leaks in the mano- 
meter tubing, whach were found and corrected during subse- 
Bent testing. 

Pamiy testmm@ata were rejected and were not reported 
Dreams Experience was galned with operation of the 
facility, the 7$ variations were narrowed to less than 3%%. 

Cylinder wall temperatures were also recorded during 
Br anne runs using a thermocouple imbedded in the wall, 
bes Econ the base. These temperatures ranged from 
Dtos” E. A plot of pope versus temperature showed 
Boseorgm latıon scatter in Py / D, ateo tenperature was less 
Since the stagnataon temperatures of the air 
aTe mere not recorded, they could not be compared with 
the cylinder wall temperatures. 

The values of Re... - 1.95x10° and Reset - 7.8x10° were 
assumed for the computed values of prae son zn rag. 20 


"Irc InodMcs he experimental values of P, / P. ISI 
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Y Zen 
X Simulated 
Heat Addition 
Š or = 2.9] 
Geh 
X 
AI 4 | | | 
HN = 1.0 To. 0.0 0,5 & 
1.0 x 
VY Computed Value 
0.9 X Experimental Value 
GQ Average of Several Values 
[me 
0.8 
0.7 
Without Heat Addition X 
x 
x X © 
0.12% x x Z 7 Xv. x Vx Vx 
= A 
> 
6 7 8 9 10 


Re,, million 


Figure 20. Comparison of Experimental aud Computed Values 
of Base Pressure Ratio vs. Reynolds Number at 
Mach 2.0 Without and With Simulated Heat Addition 
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uhiform flow nozzle without heat addition showed good agree- 
mone with the computed values. Three of the four experimen- 
tal values of Py /Po, with w/r = 3.91 compression agreed with 
the computed values. The observed value of P, / p, with 
z/r = -0.5 was less than the computed value. This differ- 
ence was attributed to an adverse pressure gradient along 
the cylinder wall that was not considered in the calcula- 
Eon In vrew Of this result, base pressure ratios computed 
by the program given in Appendix J should be considered with 
uxo che tull compression due to external burning im- 
pinges ahead of S Deseron the projectile body. 

EE Comparison Of experimental [Ref. 260] 
and computed values of base pressure ratio resulting from 
Split. 

Mimeniemorocco-Loesmneoac! there are several input param- 
R 


roms luding M, ky; k R Ee cx welatrzon- 


qu Str’ 


Pipex procedure for predicting the influence of spin on 


Cert 


P, / Po pou sussede»rewlously at the end of Section II. 
Nx Pscedure was used for the computed values in Fig. 21. 


Gieiee Ol the input parameters remains. Preliminary 


scc dSonseerescowmeuctedawath M = 2.86, RO 4 y = 1.95x10° 


and Re, e = 4.6x10°, Values of Pp/ Poo had the correct trends 
relative to experimental values. However, the calculated 
py/p, for Re, = 2.6X10° and zero spin was about 0.54. To 
decrease the calculated Dine oe the Mach number was arbitrar- 


06 
pud zero spin and" 


mi xnenoNoedetoe s. 0. For Re, 


20 D cnescalbeulated and observed vues agree, In 
, b/ Po 5 
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eO X Experimental ne 
OX Mach 2.86 
e Computed Value, Mach 3.0 
== Ey 6 T 6 
Re = 1495x310, Rennes 4. 6X10 
20 X 


Spin: 0, 20,000 RPM, 30,000RPM 
and 40,000 RPM 
0.45 


920 
30 X 
030 
Fb (40 0 
Poo 140 x 
0.4 | "VAS eeu) 
OX noe CU Ie 
30 |e 
20 X| «30 
X Nd SE 
040 ° Ü 
e 20. 20 
Ole 
30 X 
QUSS EA 40 Y 
wae 30 x 
40 x 29, 30% 0,20 X 
2:6 3.0 3.4 3.8 


Reynolds Number, million 
Figure 21. Comparison of Experimental? and Computed 


Values of Base Pressure Ratio vs. Reynolds 
Number With Spin 
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fhetrospece, additional calculations should have been con- 
meted with increased mixing, i.e. larger ku: while retain- 
Peco. Scum — 2.00. Bn tact, spin should increase Ke 
feeaisctissed in Section IT. 

A better procedure for calculating the influence of spin 
eh be to use a shift in transition Reynolds number coup - 


led with a k = k(l+(rw/U)), as discussed by Smithey, 


spin 
Nc mellcand Pubs [Ref,., 15]. The dashed segment of 
Ba! Dos vs. Re, um mns shows the influence of in- 
creased mixing (larger ky) in the vicinity of minimum P, / Poo: 

The computed values shown in Fig. 21 required four hours 
Im ceenepac IBEMIS60/67. In view of the computing time, the 
pnscedure s sugpested rn the preceding paragraph was not 
accomplished. 

In" udessemecestrimatevof the optimum location for 
Eernatoheat addition, a fixed fuel flow was MEC. to be 
mee nones ot different Length and location. The ini- 
mal adams Of the agnner streamline was 2.0 and the outer 
connie intra radius was 3.0. The heat zone were 
BESSE RK EE EE values From 1 through 6 for w/T. The zone 
with w/r = 1.0 was assigned a value of Q = eL = 0.2. 
1 c uc0359decreased by 1/(w/Y) for the other five zones. 
INEM ONES nC C positioned along the axial axis starting 
with compression at the base and moving downstream in in- 
A Sa dias.. lie results of two ol these compu- 


MOS e ern aS te two lower curves in Fig. 22. The 
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1.8 
Fuel Heating Value 
7200 Btu/lb 
ero 
Ee 
Poo 
» w/r 4,0, Q - 0.2 
w/r = 5,0, Q = 0.16 
ee 
1.0 
Dr x9 O. = 
058 
ur = 2.0 
= O e 1 
0.6 
0 1 DE Y 3 





Figure 22, Computed Values of Base Pressure Ratio at a 
Taste umocr or 2.0 Wath External Burning 
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gain in DAD for m = 0.15 lb/sec is seen to be lir] yan 
Ecuspstuve to heat zone size and placement. 

Another set of similar calculations were made based on 
Meee, o at w/r = 1.0. This value of Q corresponded to the 
= syon fuel flow rate of 0.6 pounds per second at 23,000 
fet altitude. The first three values of w/r, 1.0,2.0 or 
smüEwould noteeompute-because the.external flow was driven 
Mas onite condition by the base pressure rise and the heat 
Mao.  Ihis condition 1S analogous of thermal choking; 
however, in two dimensional flow both streamtube area change 
Ex addmwtron drive the flow toward sonic conditions. 
itemuastethree vallics of w/r = 4.0, 5.0 and 6.0, did not 
Spentenee Subsonic conditions in the external flow and 
B ursvalues of base pressure ratio are shown as the upper 
Qa C inge. 29. Ihe. high value of ney Dee = 1.91 
for the w/r = 4.0 curve was not expected. This represents 
E conso! su iicrcase im base pressure over the two 
Ewssonecurves 3n bag. 22 with a factor of four increase in 
ton Téc combustion for this flight condition, Fig. 
Boe Was assumed to start at (z/r) = -1,2 and the initial 
compression was on the shear layer, downstream of the basc, 
Po mpress tion at z/r = +2.0. There is some error 
muEcucEc s nbP3un ıntroduced by neglecting vorticity and 
Deeticmeolpless tom fan at the base, rather than the oblique 
Shock See Shaprzo [Ref. 31] for a discussion of compres- 
SAN sti aras muelles a 1" "error in the flow angle, 


Dee demps the estimate piven in Section IV, the value of 
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Figure 23. Comparison of Calculated Base Flow With And b 
Without Heat Addition at Mach 2.0 
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Pi / P remains at 1.8 and gives specific impulse of 230 
seconds. With w/r = 4.0 and z/r = 2.0, part of the com- 
pression due to combustion was applicd along the shear 
layer downstream of the critical point. With z/r = 1.0 
more of the compression was applied upstrcam of the critical 
put. but tbe base pressure was less. The variation of 
Py / P, uU /P indicates that the location of dp/ds (due 
to the external compression waves) along the shear laycr 
may be an important consideration in achieving high values 
l C Pressure wath external burning. This hypothesis 
ed tive condition where the external flow with heat addi- 
puonebecomes subsonic require additional study. 

Pure provides some details regarding the optimum 
Bor Condition shown sin big. 22 as point A. The heat zone 
era shear layer have been plotted in the lower half of 
mie co. lhe boundary layer thickness has been included 
as well as the oblique shock wave which forms when P, / P. 
IET caben Chan one. For comparison, the shcar layer 
without heat addition has been shown as a dashed curve in 
EA Pressure ratros along the shear layer have been 
Nite the upper halt of Fig. 23; again, the dashed 
e a e a ou heat addition. In both cases, with and 
without heat, the shear layer pressure ratio shows a gradual 
case distance downstream from the base. The pres- 
sıresractıo along the inner streamline of the heat zonc 
(Bommts a’ through e') has also been plotted along with the 


E NEre5PessuBe"ratros (points a through f). This 
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En camiéonesexpermrences an abrupt increase in pressure ratio 
where it intersects the oblique shock wave at c'; then, its 
Bi Sssure ratio Starts to decrease with a corresponding in- 
crease in Mach number, indicating that the flow along the 
inner streamline accelerates downstream of the oblique 
shock wave. The calculation along the shear layer was 
terminated before reaching a maximum pressure ratio due to 
Mmnemsaddie point singularity, Fig. 4, since the straddle 
method employed in the computer program had established the 
ln surred base pressure ratio within * 5$. 

Additional calculations with the computer program given 
REENEN EI J would improve the accuracy of Py, / P and would 
Meee values Lor the shear layer slightly closer to the 
ao nt but, since the critical point is a saddle 
pou eaea rity, would pot be expected to calculate through 
uS NEpornt. Calculation “through” the critical 
Petit, usime the straddle method, would require changes in 
mme program to restrict the values allowed for dF/dw, and 
force the calculation across the singularity. Additional 
aed, sicwnequired sto include this change and provide the 
s lli yo calculating the entire base flow field. 

lc x Dessscence or an oblique shock wave at 
the base, which was calculated in Fig. Du one miy consider 
the field of compression waves impinging on the shear layer 
Pvc Mmomconecentrated into a single oblique shock wave im- 
Dni on the shear layer. This simplified condition, 


oo has been studied experimentally by korst 
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feet. 3560) and was found to create an oblique shock wave at 


the base. Thus, the calculated shear layer profile and 


EN e eere ratios in Fig. 23 arc in agreement with experimen- 


tal results from the simplified condition. 
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VI. CONCLUSIONS 


LED e Od of characteristics and the Crocco-Lees 
base flow model, one can accurately predict the base pres- 
we datio Ot a two dimensional or axisymmetric body with 
moderate external burning for various Mach and Reynolds 
numbers. The er yae eaor the tyo transition 
DE molds numbers are essential for accurate base pressure 
uxo colsudatron'of a partgéscular body and should be ob- 
uuu DESEN The effect of projectile spin can 
bemineltided approximately by dividing the transition 
PoMolds numbers with one plus the spin number. 

DH uPustecan be produced by external burning. Plac- 
meethne combustion zone well forward of the base to give 
full compression on the body was found by experiment to 
paleta decrease in thrust. For maximum projectile 
Asis upersonie external flow the computed heat zone 
Mee anwas cound tO Start 1.2 base radii upstream of the 
Base and e6xtend 4.0 base radii downstream. The zone started 
Denen rsandnouter streamline of 2.0 and 3.0 base radii 
pr ecebzanmd received a fuel flow of 0.6 pounds per 
Ee Eege exhaust, ARC-168. It ap- 

p xdg rette Ste wot be possible with 
us cE D Patcsby decreasing the length of the com- 
Po crurther. Ihe £low field could not be 


Ac seat drstribution because of the resulting 
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subsonic flow and requires additional study to determine any 
wer advantage Or limits under these conditions. Experi- 
ments with actual combustion should be conducted to deter- 
mine the feasibility of the desired combustion placement 
End to measure the resulting base pressure. 

ner te zo spressume inerease along the shear layer due 
mu chesheat zome placement appeared to influence the com- 
pu devatuesof buse pressure ratio. This effect requires 
Padamtional study to determine its significance in achieving 
weer values or base pressure ratio with external burning. 
Epemnental studies are needed to determine any variation 
uw ucNcnprrrcal - - k relation when external pressure waves 


impinge on the shear layer. 
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aePENOT Aen. se DERIVATION OF THE CHARACTERISTIC EQUATIONS WITH 
HEAT ADDITION 
Equations (1) through (4) govern steady compressible flow 
eat addition [Ret. 18]. These equations neglect heat 
euet on, viscosity and body forces. In natural coordi- 
Mates Pig. t= ws. The continwity equation, Eqn. (1), 


may be written 


(21r)? owAn = constant (1A) 


where a = 0 for two dimensional planar flow and a = 1 for 
io ne Breslau. ihe momentum Equation, Eqn. (2), may be 


written for the two unit vector directions, s and n, 


Dies 
(gw xL: mLoym eod (2A) 
( ow? Ll *2P)n-0. (3A) 


iitemeierey cquation, Fan. (3) becomes 


əh 
ES DE pcs 0 
Wl (hgt > W ) = q less (4A) 


aere the total enthalpy, hg = h + aw, WaS bE cn used. The 


Be mod mame relation for entropy, S, may be written 


1 


I ae yee S Vp | (5A) 
oh 
Dre ame, Lo op 0 
Iss Ure cS) x 2 
oh 
O OW 1 op 0 | 
E (wn t p 3n? Uie DE 





Doing Eqn. (5A) into (7A) 
(8A) 


Bonanmesthe logarithmic derivative of Eqn. (1A) in the s di- 


rection 


q o «D (9A) 


Qo E 
— sin8 + — => 
T Ə s W OS on 


oje 


where the relations ər/əs = sin9 and [1/An][dAn/ds] = 90/əon 
Dave been used. Forming the logarithmic derivative of Eqn. 


(4) in the s direction 


1 3p . 130 , 1 9T 

p as p as guck ua) 
Substituting Eqns. (2A) and (10A) into (9A) 

O 1 W OW Io 00 _ 

y sing S (ce RT? s x Era oa Ü, (11A) 
From Eqn. (4A) 

q. T, jw 

T C38 REM (12A) 


MieremwEnewcenthialny has been written in terms of a constant 


EC h = Sal an vreatıne Egn. (12A) into (114) 


e ec VW ow , 90 — q G IDD a 
as RT * CST ls f 3n wc T * 077 Oe os 





Combining terms, and using the Mach angle relation cot’y = 


M^ - 1, Eqn. (13A) becomes 


el MEUM dO E usn) ` g (14A) 


76 





Thé Prandtl-Meyer Angle, v, is defined by the relation 
dy = StF aw , (15A) 


Ir nan EQ I ASE) mto (14A) and multiplying by tan u 


av 20 _ asinó _ _q 
Se tanu am tam f a an} (16A) 


Siiimerliy, Eqn. (ISA) into (SA) becomes 





oh 
IV IA 1 0 
Der as 9 Ww. on w2 ~on Su 
Adding Eqns. (16A) and (17A) 
d d E E asin6 q 
( 5; * tanu m jo 8) tanu $ T rJ 
(18A) 
f 3h. ) 
r s 
w? (“on on 
EmrsctinpoEqn. (17A) from (164A) 
Meu c n > 0) = tan asin@ Á«q 
ðS on r wc 1 
P 
(19A) 
QA. aho 
w? | ðn an] 
CO ncecharmdcreryrstic coordinate system 
DEC ES A | 20, 
secu ^q ( 53 + tank 37 DS) (20A) 
and 
ERNEUT UOS. S. E > 
Secu —3r ( ye tanu == m (21A) 


T 





With this, Eqns. (18A) and (19A) become 





(22A) 
Q(V - © : asin9 q COSU 0S oho 
- Sjny se : š > SE < 
on r s W on on 


and 
(23A) 


d(v + 0) SU sumo q . cosu 2495 | Zn 
dE T e w2 on on} 





I 


Bonus amde 25A) are the characteristic equations, with 
heat addition, governing the two dimensional planar, a = 0, 
Or axisymmetric, a = 1, flow. These equations are dimensional, 


Ceetee dimension of reciprocal length. 
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APPENDIX B: DERIVATION OF CONDITION FOR IRROTATIONALITY WITH 
HEAT ADDITION 
The equations required for this derivation are momentum 
Kaspi fo viscosity or body force, 


Dy. 1 

n ey yp (1B) 
Lucerne U 1s yelocity, p is density, and p is static pressure; 
emerey With heat addition, but without viscosity, body force, 


en zeonduetion, 


Di mop 
pare 17 (25) 


Mere ii 1S enthalpy, and q is time rate of heat addition per 
unit mass; the second law of thermodynamics 


l 


Dune E (3B) 


DN nce ysospecerric entropy; and the vector identity 


cl 


VU Vv(5U*U) - UX(VXU). (4B) 


EE bans. (1B), (3B) and (4B) 


s= * V( U-U) - Ux(vXU) - TVS - wh. (5B) 
Introducing the definition for total enthalpy, hy = h + LU-*U, 


muc cn Mbonesobtasns Croccos's Theorem 


c 


9 


TVS « Ux(vXU) - vh, * $— (6B) 
Pe eandetake the curl of Eqn. (63), 
op = VIXVS * Q-VU - QV-U (7B) 
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since V+ = V°(VXU) = 0 and the curl of a gradient is zero. 
[aie ¢b)deseripes the total rate of change of vorticity, 
qu 

Now assume that 2 = 0 prior to and during heat addition. 


Eun (7B) becomes 


= = VTXVS = 0. (8B) 


ImenssEanz (8B) is the condition which must be satisfied dur- 

a C on the flow is to remain irrotational; 

mo e temperature pradient must be parallel to the 

entropy gradient. Referring to Eqn. (3B) and assuming 

h = Seb where Co Bes came, che condition for irrotation- 

ality may be stated in terms of pressure; that is, the tem- 

puc NPEddmNent must be parallel to the pressure gradient. 
oo steady flow, the irrotationality condi- 

Aa Oe eien in terms of heat per unit mass which has 

been added, Qaz In Eqn. Ceea = UO, and, using the 


total enthalpy h, = h + LU-U, 


EC - Q) = 0. (9B) 


Menal ong any Streamline 


ho =O" = In, (10B) 
where it is assumed that the upstream conditions, before heat 
SO are uniform (Ret. 19]. With these as- 
sumptions one may take the gradient osram. (10B) 


vie - Q) = 0 (11B) 


§ 0 






since ho sin eche steady form of Eqn. (6B), 


00 


Teche _ 0, and Eqn. (11B) 


VS = 


Vh, - 4 VQ. (12B) 
Deo, from Egn. (11B), one may write 

VQ = V(h + 4U-U) = e, VT EEUU | (ESB) 
emong Eqns. ('Z2B) and (I3B) in (8B) 


e [va - vm) X T VQ = 0. (14B) 
D 


nee (VO)X(VO) = 0, Eqn. (14B) reduces to 
(VQ)XV(5U*U) = 0. (15B) 


Par. tnC condition which must be satisfied to pre- 
co tationality 1n steady flow with heat addition; that 
meee Orddment OQreeneat added must be parallel to the gra- 
enim kinetic energy. 

Tr the gradient of heat added is a small perturbation, as 
Rms mGase Omelincarrzed flow, then the irrotationality con- 
Bao 5b) sesutomaticaliy satisfied. To show this 
U' 
U 


oo 


M 0 - hn20', where pu and 0 = U (i + We where g'eg! 


EPUM ubstepcyuting 3n Egn. (15B), 
uos a Q Y(U yUt* + %U-U') (16B) 


NNUS coded third onder perturbation. 
This Condition may be applied to Eqns. (24A) and (25A), 


Appendix A, and neeslect the vorticity term 
1 [38 | aho 
W on on Í 
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o 
Pig samplityime the method of characteristics calculation. 
Buc SII. the heat zone is started in a region of 
Mito rato along a line normal to the streamlines, and if 
the heat is added uniformly in a zone bounded by two stream- 
Buchen the Condition of Eqn. (15B) will be approximately 
Satisfied within the zone since the gradient of heat added 
will be approximately parallel to the gradient of kinetic 
cos Second if the first condition is not satisfied 
Me, alene the streamline boundaries of the heat zone where 
Mieremts a jump in QO normal to the streamline), the vorticity 
ms tillbe neglected 1£ the heat added is a small perturba- 


tion, that is 
q = U-VQ = h,U-vaQ'. (17B) 


Then the heat term used in Eqns. (24A) and (25A) must be 


h  w(9Q'/3 
eee Ios Eo 1/r,. (18B) 


wc. T wc. T dS 
p p 
MNAE ion 15 Satisfied for the heat zone boundary, where 


cr EN INI UU to pet - 45 deflection of the 


Mier heat zone Streamline. 
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APPENDIX C: BRIEF DEVELOPMENT OF CROCCO-LEES BASE FLOW FORMU- 
LATION [Refs. 1,8,13] 


The mass flux, m, in the base flow is given by the inte- 


gral 


6 
m A (2rr)” pudr (C) 
0 


where 6 is the viscous flow thickness, r is radius, p is den- 
wees axtalecomponent of velocity, we equals zero for two 
dimensional planar flow and «a equals one for axisymmetric 


flow. The momentum flux is 


I = ^ (rr) pu?dr . Ceo 
0 


For the base flow, where there is no wall drag, the momentum 


equation becomes 


dle- dm _ D ClO. 
Hu m (T$) 6 SE (2C) 


where u, is the external axial component of velocity, dp/dx 


Bee Te pressure syadıent and, assuming small angles, 


m dé 
nu (276) Bee): (4C) 


e 
oia the Grocco-Lees theory is that 


dó 

2a = 5C 

dX a k (5C) 
Wiete sh is a ton dimensional mixing coefficient which must be 
speciited. the value for k is approximately 0.0033 for laminar 


Dt for turbulent mixing. 
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"A mean velocity is defined by 


u, = I/m (6C) 


wea Velocity ratio by 


K = u,/u,. (7C) 
The mean density, 01» js sesio Utic equatuon 
m = (16) 6p,u, . (8C) 


assumed that the stagnation temperature is constant 


pm os Oe i On and that the lateral pressure gradient, 


r r iS ero. Knowing the pressure, p = Po? and Oy» A 
mean temperature, T4, can be defined 
T, = po / p, R. : (9C) 


Ian incinse the displacement thickness, ó*, and momentum 


Bmrekness, O°, On ordinary boundary layer theory 





ó (10C) 
re). FF enn. D jar 5040 m 
0 Pole pu Tre 
e e 
and 
6 
EOS Io S. E _ u 
DE - SL Gu Gees (1 a 
0 eye e 
(11C) 
¿(1+0) I T1020) _ ] | 
o_u?r” 
e e 


Notang ihat k = I/mu, sur Sbanse (10C) and (11C), one 


Das 
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s ta) = gx (1+0) Y sss (1*0) 


K = : Fa Ben IE . (12C) 


BEE Et transformation [Refs. 8,13] 
e os soc, aideso-* in Eqn. (12C) are replaced by the 
incompressible quantities ò., $: and 627. Thus, within the 
Fes Urietilons ol the Stewartson transformation, any compres- 
sible laminar base flow can be related to a corresponding 

M conmpresswble laminar flow. It ais assumed this also applies 
in the turbulent case. The ratio of mean temperature, Ti» 


Moms baciannon temperature, 1% may be wrıtten 


JE 
ME et en 
T fe 5 Wy (MSc) 
S 
where 
y sien) 
f = 14C 
sta) . gs (1*0) pes 
i i 
DES ICH and a. estra enation specd of sound. It has 


K E imo rc conventent to work with the quantity F > 0, 


E + 1 = f/k?. (15C) 


Experimental data [Ref. 8] has shown a relationship between 


F and « which can be expressed by the polynomial 
Ore, tor - 0-0T185F*, (16C) 


IS ic 0 72-amd dk/dF = 0.0. 
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From Eqn. (8C) one may write an expression for mass flux 


a 
u Y(T6) 6 pue 
m = ma_ = — (EC) 


; mer 2 
AM TO! 5 Wo] 


Usame this and Eqns. (3C), (4C) and (SC) one may form a 


Bierce Ordiliary differential cquation 








A len 
z E 3- WE I BEC! E AE 
W - -> 
e 
| NE (TSC) 
1 —— 0D 
KF 2 e e 
e f DN ND Pos ay k 
1 





2 
° 

s 6 

w. (1-Y 12) Iu SG gë mis (de /dF) - ? 

Eqn. (18C) has a saddle point singularity, [Refs. 1,8], also 
uL e Critical point, which, in integrated form, detcr- 
mines the one unique solution describing the given base flow 
condition. In two dimensional planar flow without external 
heat addition, where 9. is known from pe pt ne Prandtl- 

MEE aio one ma start at the critical point, where 
Bechzchemumerator and denominator of Eqn. (18C) are zero, 

Ai teo pated toward the base [Ref. 8]. In flow with ex- 
ternal heat addition, and in axisymmetric flow, D. is not 
known as a function of We thus, in these cases, one must 
ao Eon. (18C) at the base and concurrently 
Sonve sehe 2 tersıtal flow t1erd, by the method of character- 
ileso Fımearized theory, for the values of D. along the 
Armee iiss process 15 one ot trial and error described 


i section iV as the straddle method. 
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APPENDIX D: CALCULATION OF BOUNDARY LAYER MASS FLUX AT THE 
BASE 
The boundary layer is assumed to be adiabatic, have con- 
EE across its thickmess, ó, and have a velocity 
profile (u/u,) = Bon nu where U, ie mete rial i low 
May and n s an integer [Ref. 24]; 1 <n < 7, for 


fainter eto tunbulemestlow. The boundary layer mass flux at 


sbs 1S Given Dy the integral 


E FO 
E b 
my, = J (271) “pudr (1D) 
b 


where a is zero for two dimensional planar flow, or one for 
axisymmetric flow, p is density and zs reece radius ©, 
Epsum sine the equation of state, p = pRT, and the 


aove assumptions, Eqn. (1D) becomes 





Ee a Zr 15m 
p (206) usv ) dr 
c C ee S 
IDE T eee G ss s s (2D) 
b R u2 DH 
d P e ( by 2/n 
b T 2 Š 
p 
where Ba "NM Een pressure. Defining Wo as u/a., 


where a. ENEE peor of Sound, Eqn. (2D) becomes 


18 


a a b 
Dn M de Y6p,w, (276) $ (— 


. Š ) 
£ 
b el 
DE) 


MES (3D) 


 — — — —— 


LEST ` 
wê Co 2/n 
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Then 


1 r 
m, E O^ ax 
ma E E ey 
yop w, (216) 1 - B?X 
0 


= = Zah > 
where x = (r r,)/6 and B 5 Woe 


If a = 0 and x = SE Eqn. (4D) becomes 
Í 
n 
_ zmidz 
eer of eae sr 
0 


I emo be integrated. Form odd, Eqn. (SD) gives 


(0E 19072 
1 1 d 
I = n 4 ——— ln(——) - p 
Dopp a ne^ E [n+1-2j]B J 
F e umma lon 1S omitted if n = 1. For n even, Eqn. 
(SD) gives 
n/2 C7) 
1 TEB il 
I = n 4 —— ln(— ) - A ee 
2DEVN par 1-B A (nei - 21] B2? 
If a = 1 and x = m Eqn. (4D) becomes 
na Ty Do n, 
Ero f 1-B?;: " 6" $ 1-B?z? e 
0 0 


I S i r JM Bano (8D) is the same as Eqn. (5D), 
ie S i T C Cin Equs. (6D) and (7D). Integrating the 


first integral, Eqn. (8D) becomes 
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n 
` 1 CB 
Lut "P Eus In (7) 2. 


a EP 


r (9D) 
TM 
DNO D. 
Mac sesuibseosmiens. (6D), (7D) and (9D) are used with Eqn. 
(4D) to evaluate m, or m, . This boundary layer mass 
2D AXI 
flux at the base, before the corner expansion, My > is set 
Simao cre mass flux after the corner expansion, 


Mp > Ende 


ENEE temes to the Crocco-Lees formulation, Appendix C. 


6 I ve Eet (10D) 

m = On W = m! = —rnM—o V 

SEH m Way e! [F+1-L hy 1 2] 

remo ge the value of w. assumed from the straddle 


Doea one Solves Egn. (10D) for the value of F' which starts 
theontesration along the shear layer. Similarly, in the 
axisymmetric case, 

yo ew 

8 nt p _ ee 

my, = 2TYŠ PoWolayt my ————————— (ITD) 


AXI Gee s Lat 
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APPENDIX E: MASS FLUX CALCULATION ALONG THE SITEAR LAYER 


PaO sendet tan, (17€), the mass flux, m, may be 


written " 
"ol 8 p,w, 
a ee (1E) 
K[F+1-5 w2] 





where a = 0 for the two dimensional planar flow and o = 1 for 
Gu imetric flow: One may evaluate Eqn. (1E) at two adjacent 
Points along the shear layer, m. and D. 1 and, solve dor the 


same between these two points, dx, by using the additional 


a lacions (Ref. 1,8]. 


dé 


and 
dm . Ol 2 dó ` 3 
no (To) PW, a. ( dx Quy (3E) 


Aene aenal l anelles, Eqn. (2E) states that the rate of 
eremeenon tie Shear layer height or radius, dé/dx, is equal 
O e en anelee O plus some constant, k. Then the rate of 
GREEN EECH by Eqn. (3E), is the mass 
cine ele shsar layer from the external flow. 

To compute mi,» one increments F, with AF, , Say AF. - 


0.1F, with the genden Lon AF 220,08, 


E = + AP... (4E) 


W iy m AY; / (dF/ dw.) (SE) 
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where dV / dw, Pmconpipeaetnomecan. (18C). The value for SES 


is computed from an empirical relation using Fix given by 


bane (16C), with Ki IZ a F. 5. Wee 1Ssentropic 


+] jsp 


relation between static and total pressure, in terms of y 

and Wo is used to evaluate p, at i and i*l. Average values 
between i and itl are used to evaluate dé/dx and dm/dx, where 
bic 1sentropic relation between static and total density must 


memused with the latter. Since 


m. 
1+1 


m. + (dm/dx) dx (4E) 
and. 


BC? 


6. * (dé/dx)dx . (5E) 


Ma solve directly for the value of dx, in the two dimen- 


Sional planar case, which corresponds to the increment Ar; 





m. m. 
pr dó dm 
dx, = ó. (— - #— )/( S“ - S). (6E) 
Lon ] dÉ 8:411 dx dx 
For the axisymmetric case 
1 dc 2 
dx. = > ( dx. - — += dx: Je (7E) 
ES —— S. ar 


Pan m UE may be iterated several times, as long as dx. 
AXI 
IE oE I pimi ne the correct value of dx. 

` E 
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Serer, Fs  CALCUBATION OF BOUNDARY. LAYER THICKNESS AND 
ME ENGT COEFFICIENT | 


The laminar boundary layer thickness, ó was computed 


pe 


m 5 J n9he Compressible Blasius solution [Ref. 8] 


6, = (5.0 + 1.2 (y- 1)M*) x/ (Re) * (1E) 


punc ae ance along the surface from the leading 
edge and Re, 1s the Reynolds number based on this distance. 
If the base Reynolds number, Re, , is less than the transi- 
Crona Reynolds number, Re,» Licimerce boundary layer at the 
base is completely laminar and its thickness was computed 
icon with x taken as the body length. 

If the base Reynolds number, Re, , Peer cai r than the 
ani ron Reynolds number, No the boundary layer will 
De partly or completely turbulent. The turbulent thickness, 
On, was computed from [Ref. 8] 

5 3e y 178 Cen! Cea) 6/88") d 
Cei Nu axmanusemgmconpressrble turbulent skin friction 


coetticient given by [Ref. 24] 
4 > 
(0.242)/ (Ce; ) = logy, (Re, Cy; ) (EF) 


and Com is van Driest's compressible turbulent skin friction 


cocus xuEcrvencby [Ref. 24] 
| (4F) 
ue s " 2 
eB Sloe, (he Cag) + 1.261op, (1-B2) 
(BC y 10 x fM 10 
fM 
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where p^ 1s (v-1)w2/2. Egustrons (5P) and (4F) must be 
solved by iteration. The quantity 6/6** (Appendix C) may 
Da aii eeecadmr ore tions two dimensional 1/7th profile using 


the same assumptions employed in Appendix D. Briefly 








b 
k = (pu) ts 
T y as (^d " )dr (SP) 
b À NDA mo (6F) 
Ber E ule N ( 5 ) B ( Š ) 
Š Z e d CELE (dr/6) 
b EC <6 





-= = ? (1- B^) (4 m - l 73 
2B 1-B i 3(8-2j)B*- 
j=1 
A 
s gn T d —— ). (7F) 
2B iz (9-23)B 


ENIM nce cuwmalMi5insrsftion process is more or léss probabil- 
:stic in Nature, it was assumed that the transition from 


laminar to turbulent flow can be described by the integrated 


normal distribution. It was necessary to assume a turbulent 
transition Reynolds number, Re rp? formers purpose, where 
Re, + Peep r oximately four times Ko, The transition 


Gum ias made with the error function, erf, which is 
Ente Toschesinteprrated normal distribution. The 
arrument of the error function was formed in terms of base 


incu cneehesinPtralstransition point, and was weighted 


d 





so its value was -0.96 at Re. and 0.96 at Re ert: A numeri- 
cal span, S, of 4.0 takes the integrated normal distribution 
TEENS MIOSUMSSPormbhatrng these conditions in terms of 


body radius, r, 


A Re O) (Up) (8F) 


K I ans ition point on the body measured from the base, 


and 
Lert = (Re... /Re, ZO ODIT) (9F) 
Cea Duon transition point. Note that arr will lie 
dowmstream of the base if Re, t > Rey Forming the average, 
ceca oe ¿ eo (10F) 


ee Lance from the base to the center of the transi- 
Alone one £ e NEE EE tire “cumulative 
a O Ust be 05. The value of erf at this point 

O and e cumulative distribution integral, Icp> 
is 


I 2 wu sts emp. (11F) 


CD 
The values of erf at the base, ERF, > range from 1.0 to 


O Ae rms OZ 





ERE, = erf(— TE (12F) 


Ihen, when Re) > Re, ÖBASE 1s given by 


á = 46, (1 + ERF,) + ro (1 - ERF 


BASE p? : 
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intent lor), ERE will be positive if Re, is less than 


b 


thus as Re, ranges from Re,.. toe thc 


L 


C 


AMINO LI range from 0.99 to 0.01 and the con- 


L 
rau oT od d e Om 0.01 to 0.99. Pie valwc for 


(meri to compute the velocity profile at the base (see Ap- 


pendix lejeers formed exactly as Eqn. (13F) 


= 4(1 + ERF,) + 4 (1 - ERP (14F) 


DBASE p? 


mince n ranges trom l to 7. 

When the transition zone lics along the shear layer, 
GEHEIERT a transition value for the mixing coeffi- 
cm wk cat the point being considered, Xpj Ihe distance 


X is measured from the base in terms of r. To carry out 


BL 


wwe calculation, the value of erf at Zei: PRP oi: EIERE 


EE e 1n terms of Z 


S UC x CU 
m MEE . trav BI 


BL 
2 Tort I Zuel 


(15F) 


te 


Natera an (15F) réduces to (12F) since Xp], 1o zero at 


EE Eeer Gresch the valuc for the coefficient, 
Ko: JS oI cina imi lar manner trom the laminar and 


bulent mixing Coefficients, k, and K... 


L 


+ ERE.) +aDKk (l - ERF (16F) 


BL BL? 


Jusi aa the value of ERE will be positive if Z 


B I: 


aya Vice. Versa. 


trav 


1s downstream of Xp» 


GES 





APPENDIX G: DERIVATION OF LINEARIZED STEADY FLOW EQUATION 
WITH HEAT ADDITION 


ux ccuscrons 1 (4J). are required for this deriva- 
tlon 

Ve (pÚ) = O (1G) 

pU*VU = - Vp (2G) 

U-V(h + 5U*U) = q (36) 

p = pRT (4G) 


Forming the directional derivative, U*V, of Eqn. (4G) 
U-Vp = pRU-VT * RTU-Vp. Di 


Eet Ean. (1G) into (56) 


pRU-VT = U-Vp + pRTV-U (6G) 


ENEE of Eqn. (2G) in (6G) 


RU-VT = - U-v(4U-U) + RTV-U (7G) 


Assuming = is constant, using h = Se and SE = YR/(Y-1), 


substitute Y/(y-1) times Eqn. (7G) into (3G) 


1 


VeU - a UN UU) -= (Gef: (8G) 
a 


INN ONU EU Uh. WsUuY-ve. 


With this, Eqn. (SG) becomes 
(9G) 


(URT Vp). V 
(4U2 o Hä, nr 5V$-V6) 7 q'/(1 * T!) 


DUE 
ECL T!) 
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[eni oon] y irst Order terms, Eqn. (9G) becomes 


2 e 2 = t 
UE a. (10G) 


This is the linearized steady flow equation, with heat addi- 


ron Fox sırrotatıon subsonle or supersonic flow. 
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APPENDIX H: DETERMINATION OF THE MACH ANGLE 


Inti derives a numerical procedure for deter- 
mining the Mach angle, u, for a given value of the Prandtl- 
Heyer angle, v 


e aa E Meyer angle may be written [Ref. 18] 


l m2 -1)? (1H) 


c 
tl 


- 1 = 
° tan lo (M2-1) 5) - tan 


mere M is the Mach Number, 1 < M < e, 


((¥-1)/ (y+) )7 (2H) 


m 
II 


TOR iS the specific heat ratio. Let 


uot € = 1/2 (3H) 
where u is the Mach angle defined by 
sin u = 1/M. (4H) 


Dom Eqns. (3H) and (AH) 


(M2-1)% = tane. (5H) 
Wenn requatiuon (SHA), equation (11) may be written 

tang(vte) = gtane. (6H) 
im cio relation between v and € is still implicit. 
For y = 5/3, g = and Eqn. (6H) can be reduced to the ex- 


plicit tomm 


tan 


Nim 


= tan; ( 711) 


i" oA oj) and (71), the simple relation may be 


written when y = 5/3 
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IER 2tan ! ((tan7) ). (8H) 


iimcwoialyatcal result (Ref. 37] was first given in terms 

of M, rather than u, along with a more complicated expres- 

A A There are no analytical expressions, 

u = u(v, y), for other values of Y, consequently, u must be 

atea mincd amen tables, ser numerically, in all other cases. 
Equation (6H) may be written in a form for iteration 


which always converges to the correct value of £ 


B -1,1 
£:,1 = tan PE (9H) 
o f (p c-r 2 and 0 <v < y . The value v is given 
- i - - = max max 
my 5 (Ti) tor M infinite 
TT L 
"max ^ S m - 1). (10H) 


One may start the iteration of Eqn. (9H) with an approxima- 


tion for £1 gv n boy Eon. ER 


S 


e 2tan ! ((tanY) E CET) 


This approximation fails for large values of v since it may 
give a value £1 Poy 2. einen this occurs, set Ey * T/2 to 
Tohe iteration with Eqn. (9H). 

LINIEN III always converge to the correct value 
omer meaivereence 1s Slow and, after two or three itera- 
ps aula be replaced by a method utilizing the derivatives. 


Ime ance Hs used to find an approximation of e 


Uche y se Of Convergence of a faster method. 
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One method which provides rapid convergence may be written 


ONO by using a Taylor's series expansion for the 


separate sides about Es and solving for Ae. Let 
f - tang(v*ei) (12H) 
n 
f, - gtane; (1511) 
1 
then 
Ag N 
SEI =) a + o e: cort 
1 l 1 1 
REN) 


Euer ER - tr) 
1 1 d 1 


where Ac, = Uris he tterattonm in Eqn. (14H). The advan- 
Pore OL mms form lies in the fact that the derivatives may 


c. Tt! 
1 dn E 
and fj involve the See xs whieh em Demwritten as 1 * tan*x, 
1 


n the derivatives may be included with only two 


Bee written in terms of the values of fi and f 


noc e va wat Tons by the computer of the tangent functions 
A Eeler Depending on the'number of deri- 
Moves included two Or three iterations with Eqn. (14H) 


Mould De wstrticlent. Them of course, 
p wc o - NES (ARS) 


is used to evaluate f, and f, E 
1 Je] 


TSN KEES the above procedure will 


Ges ony = 1.4 and 


provide twelve place accuracy for g in two iterations. With 
the required value for €, uy can be written using Eqn. (3l). 
This method was used in the computer program discussed in 


Appendix J. 
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APPENDIX T: SPINNING CYLINDER DESIGN 


The spinning cylinder was designed to be compatible with 
n o i li co jet fixtures, nozzles and spacers used for 
I -EncuspaDnLINPEUDVestipatu)on [Ref. 29]. The spinning cyl- 
imder and its support assembly had the same external dimen- 
Sions as the stationary cylinder design, 1.282-inch diameter 
wa S inch length, to insure similar test results without 
“pene The cylinder was driven rn ini ral impulse tur- 
bine. The design speed of 48,000 revolutions per minute, to 
naccherhe "span number ot 071595, was not achieved. This was 
due in part to physical interference during the fabrication 
Process wiich reduced the tangent weds component of the tmmbine 
Heuston lessetian one-half the desired value. The maximum 
speed achieved, without nozzle flow, was 18,000 revolutions 
per minutek 

The spinning cylinder was made of Plexiglas, to reduce 
Pottoman was glued to a metal endcap which formed 
tho inden sna em incmciidcap was driven through a shaft 
which screwed into the internal turbine. The shaft also 
EE EE n anal alr bearing surface 
OI EES with pressure and friction 
tonc e nozzle flow. "Pis part of the design was 
dee e e a the leads due to nozzle flow caused the 
tud Elend stopped the cylinder. Time did not 


Pica tacon CO correct this deficiency. 


AL 





The spinning cylinder was positioned radially by 24 air 
bearing pads provided along the outer surface of the support 
cylinder. High pressure air was supplied to these pads 
through orifices formed by drilling into supply tubes sol- 
dered into the walls ot the support cylinder, Fig. 11. The 
bearing air exhausted through ports to the inside of the 
support cylinder where it mixed with turbine exhaust and 
flowed out the bottom. 

The base pressure on the spinning cylinder was measured 
by a total pressure probe, positioned on centerline next to 
the base, and a mercury manometer. The cylinder speed was 
determined through a photo diode and a pulse counter. The 
photo diode was illuminated by a laser beam reflecting from 
the metal endcap and half of the endcap was painted black to 


cenerate one pulse per revolution. 





Disassembled View of Spinning Cylinder. 
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APENDIX J: COMPUTER PROGRAM 


The FORTRAN computer program, included at the end of 
is p ri required I50K for compile, 110K for link and 
Actor execution" If was written to compute base pressure 
NEO. within approximately 0.5 percent, by applying the 
equations and methods discussed in this thesis. In view of 
Pree cCiipirical information and approximations employed in 
many of the equations, computation to greater accuracy did 
not appear warranted. Using an IBM 360 computer, computation 
of the base pressure ratio for one flight condition required 
Epox mately two Minutes for the fully turbulent shear 
meer Without heat addition. The addition of heat upstream 
of the base and low Reynolds number, which involved a long 
Misc he ar bayer, increased the size of the characteristic 
mesh and resulted in computation times as high as ten minutes. 

The flight condition was specified before making a com- 
Mas aneluded choosing the Mach number, card 80, 
and the base Reynolds number, card 140. The two no spin 
transition Reynolds numbers, cards 120 and 130, were choosen 
Immune uEDar body to provide agreement with available 
I jr n i I he tratio of body radius to body length 
was specified on card 150. The spin number, rw/U, was given 
E wd ine Tarıo of speeitie heats for air was set at 
IN Cara too and 190 were interchanged for two 


A ana low) For flow wathout heat addition, 
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AUN as set to Zero, card 60 to ISTRT = 90, and card 100 
BONS TRI = 20.0. 

Several additional considerations were required to com- 
Dor eee ente condition with external burning. The start- 
tes Tor the radii of the heat zone inner and outer 
streamlines were specified on cards 420 and 430. The axial 
positions to start and stop heat addition were given on 
En DUESMdCVU with reference to the base. The value 
ROMO mmcCards 90m was specitied to obtain a desired - 4%” de- 
fiection of the inner heat zone streamline, however, this 
value could have been specified by assuming various fuel 
flow rates and heating values. Finally, the value for ISTRT, 
card 60, wac set ' low enourh (minimum value of 1) to place 
the Mach line boundary condition well upstream of the outer 
Sa ee on waves generated by the heat zone. The values of 
EE deduantitues were printed, to monitor the progress of 
Ene computation, in addition to the required values for the 
base pressure ratio. 

A simplified flow diagram is shown in Fig. 1J. It is 
ap endecduro-shon the basic sections of the computer program 
and computational details have been omitted. 

In the computer program the group of IF STATEMENTS which 
ice sh pointe for heat addition can introduce a 
in iss E I two corners of the heat zone. The amount of 
Mane teen two points is proportional to thë ele- 
DUDEN solos the eharacteristic in the heat zone. 


(ieee SE in radial position against the 
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EEES ETCH CONDITION VARIABLES: 
M, Re,, Re, Re, i» Y, ku, Kr, pg. Heat 
io es mm E uL start tor 
Characteristic Mesh Upstream of Base to 


Enclose Heat Zone (asp eee = -0.]1(Jk00-1)). 


Gompute-Base»Boundary Layer:  r/ó and n. 


Ponpuces characteristic Mesh Along Body and 
Ie zone Up to the Base. Store Values 
Along Base Characteristic to Reset for the 
Straddle Method. 


Compute Mess Flux at the Base. 





oj Gn Value of Wo lo (CT ıcerlles Netlog, 


Compte Expansion Fan and F. 


Conmpurteminmcrenent Along Shear Layer from 
der p. dw. = dF/ (dF/dw_) and 


mur Du" (dm dax)dx. 





Compuce Characteristic from Shear Layer to 
Up ar malacehelkıne. Correct @ at Shear Layer 


Poceo e one Shear Layer Until Critical Point 
Where dF/dw, Becomes Negative. Compute pA 


from Value of wg at Base. Assign Next Value of 


w. at base from Sign of dF/dw, Numerator. 





STRADDLE LOOPs 


STORY 
Eu xU ili ied Flow Diagram for External Burning Assisted 


Projectile (BAP) Computer Program. 
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‘Inner and outer streamlines of the heat zone, and in axial 

po metonmacatmist the stirt and end of the heat zone. Thus 
element rat the lower left or upper right 
€COimer OL the Meat zone could be outside the zone, but 

still pass the IF condition tests and be assigned % value 

ETE cam es error is small with a small characteristic 
Ian cis selected. "The same applies to a left running 
Since upper ett or lower right corner of the heat 


ZONE. 
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